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GENERALIZATION AND REFINEMENT OF AN AUTOMATIC LANDING SYSTEM 


CAPABLE OF CURVED TRAJECTORIES 

Windsor L. Sherman 
Langley Research Center 


SUMMARY 

Refinements in the lateral and longitudinal guidance for an automatic 
landing system capable of curved trajectories have been studied. Wing 
flaps or drag flaps (speed brakes) were found to provide faster and more 
precise speed control than autothrottles. In the case of the lateral con- 
trol, it is shown that the use of the integral of the roll error in the roll 
command over the first 30 to 40 seconds of flight reduces the sensitivity of 
the lateral guidance to the gain on the azimuth-guidance-angle error in the 
roll command. Also, changes to the guidance algorithm are given that permit 
•TT-radian approaches and constrain the airplane to fly in a specified plane 
defined by the position of the airplane at the start of letdown and the 
flare point. 


INTRODUCTION 

An automatic landing system capable of guiding an airplane to a precise 
landing would be of distinct benefit to both civil and military aircraft 
operations, as the cessation of operations due to inclement weather would 
be greatly reduced, as would the use of alternate airfields, which are 
often a long distance from the original destination. Reference 1 describes 
and reports the oerformance in a bland environment of an automatic landing 
system that can guide an airplane along a steep, curved trajectory to a 
precise landing without the use of a nominal trajectory. Reference 2 con- 
siders the same system in a nonbland environment. The disturbances consid- 
ered were wind, wind shears, turbulence, data sample rate, and control- 
sur face-actuator natural frequencies. Data rate compatibility with the 
microwave landing system was shown. 

The work reported herein generalizes the basic system described in 
references 1 and 2 so that the guidance commands can be generated from 
data acquired by ground-based acquisition systems. Constraints on the 
flight-path angle and the use of aerodynamic surfaces instead of automatic 
throttles for speed control are investigated. Modifications to the lateral 
guidance system that remove the sensitivity of lateral guidance to the 
gain on the runway centering error are given. Unless otherwise noted, 
all data presented herein were obtained with a steady wind and the patchy 
turbulence described in reference 2. 


SYMBOLS 


The International System of Units is used throughout this paper. 

All angles are measured in radians. 

D operator , — 

dt 

g acceleration due to gravity 

h altitude 

kj G integration gain in lateral guidance 

. . .k^g gains 
m mass 

q pitch rate 

R range, slant distance from ground data station to airplane, posi- 

tive from airplane to radar 

R h projection of R in horizontal plane 

R-j >| component of parallel to runway center line 

R -12 component of R^ perpendicular to runway center line 

R ^ 2 component of R in vertical plane 

T thrust 

t time 

u acceleration in x body direction 

U w ’ V w ’ W w windspeeds along x- , y-, and z-axes, respectively 

V Q speed command 

V T speed in inertial space 

w vertical touchdown speed 

z j 

X,Y,Z inertial axes 

x a R^ - 900 . 0 

x ,y , z coordinates 
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I 


a 

Y 

6 f 
6 SB 
e f 
? 1 

C 2 

e 

6 a 

0 AA 
9 AC 
9 e 
5 


* 


angle of attack 
flight-path angle 
flap deflection 
speed-brake deflection 
flap angle change 

runway coordinate, positive from landing point to opposite end 
of runway 

runway coordinate, positive to right of runway 

pitch attitude 

radar azimuth angle 

actual azimuth guidance angle 

command azimuth guidance angle 

radar elevation angle 


= tan 


-1 


h - 20 


l 11 


- 1175 


servo time constant 
roll angle 
heading angle 


Subscripts : 
c command 

i initial 

o output 

p plane 

A dot over a variable indicates differentiation with respect to time. 


BASIC AUTOLAND SYSTEM 

The autoland system used in this study is the system described in refer- 
ence 1, modified for the effect of wind as described in reference 2. In 
the system described in references 1 and 2, an airborne radar was assumed 
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and modeled as part of the system. Actually, the only data needed are the 
components of the range and range-rate vectors along and perpendicular to 
the runway center line . It was assumed that the airplane would be supplied 
with that data. Several methods of using these data in an airborne processor 
to obtain the inputs to the guidance laws are: 

(1) Use range data to obtain \|>, e AC , 9^ , and y q . 

(2) Use range data to obtain ®aa ’ anci °btain Y Q and ’t' from 

an inertial measuring unit. 

(3) Use the range and range-rate data to obtain ip , 6 ^q, 6 aa> and Y c • 

The first method was discarded because small differences of large num- 
bers were involved in the determination of 4>. The second method is essen- 
tially the one used in reference 2. The third method was used to obtain the 
data for the guidance laws in this paper. It was assumed that the three com- 
ponents of range and range rate , R-] 2 > ^13 anc ^ ^11’ R i 2 > anc * R 13> 

respectively, were supplied to the airplane guidance system. Tne equations 
used to determine the information for the guidance laws were 


R h = 


+ r 12 2 


(1) 

R = 

^7~ 

+ R>I2^ + 

Rl 3 2 

(2) 

% = 


+ r 12 2 


(3) 

R = 

K 2 

+ ^12^ + 

C\J 

on 

(4) 

0 A = 

: sin - ” 1 

( R 12/ R h) 


(5) 

6 e = 

; tan - ”* 

( R 13/ R h) 


(6) 


(See fig. 1.) The heading angle of the airplane was calculated by using an 
inverse tangent method based on the sine and cosine so that the correct 


quadrant was obtained. The equations used were 

sin i|) = R 12 /$ h (7) 

cos i|» = R 11y /ft H (8) 

and the heading angle was obtained from 

i|) = tan“”'(sin lp/cos i|0 (9) 

This is equivalent to determining ^ from tan -- * 1) because the 

inverse tangent is defined only between 0 and ±it/ 2. The use of equa- 
tions (7) to (9) removes ambiguity as to the proper quadrant for the angle 
4). The other data were calculated as follows 

e M « “A - * (10) 
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( 11 ) 


6 AC = tan_1 


100.0 

COS 0 


AA 


If the initial heading angle 4^ is greater than V2 , the angle 
passes through tt/ 2 during the landing maneuver. This causes the tangent 
°f 0 AC to become infinite and large roll angles are commanded. Therefore, 
if 4^ was greater than ^/ 2 , 0 AC was calculated by the alternate 
expression 


6 AC = 02) 

throughout the landing maneuver. In addition to the angles 6 .. and 640 , 
the commanded letdown y ’ and the flight-path angle of the airplane y 
were needed for the vertical guidance. These were given by 


Y c = 


-tan 


-1 




R n - 20 


11 


- 11 


75) 


2 R 2 
+ r 12 


Y c = 


• -tan 




R n |tan y 20 | 

20.0 


(H > 20.0) 
(H < 20.0) 


(13) 

(14) 


where Y 2 q is the value of the airplane flight-path angle at an altitude of 
20 meters. The flight-path angle of the airplane was determined by the use 
of the expression 


Y 



(15) 


The quantities 0 ^q , 0 ^^ , and \|) given by equations ( 9 ) to (12) were used 

in the turn algorithm, which is 


*0 = * 2 [*e - * - k 3 » * k, (8 #c - e #A ) 


( 16 ) 


where is the heading of the runway and the gains and kp were 

calculated by the method given in appendix B of reference 2 and is 

a constant gain and is equal to 0.8. The quantities y c and y, given by 
equations (13), (14), and (15), were used in the letdown algorithm to gener- 
ate a pitch-rate command for the autopilot. The algorithm, given in refer- 
ence 1 , is 



g tan <l> 


sin 


+ a 



(17) 


The gains and conditions on the use of the letdown algorithm are given in 
appendix E of reference 1 . 


Winds and Turbulence Used in Study 

Steady wind and turbulence were used in most cases presented in this 
paper. In the program used to generate the wind and turbulence, headwinds 
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were negative and down blasts were positive. The steady windspeed was 
-25.8 m/sec and was directed at 77/4 radians to the runway center line so 
there was a component of -18.24 m/sec along the runway and -18.24 id/ sec 
across the runway. There was no vertical component of the steady wind. The 
turbulence was based on the output of a random number generator that gave a 
normal distribution of random numbers between ±1 . The specific turbulence 
used was the patchy turbulence used in reference 2. The turbulence was cal- 
culated along all three inertial axes, then summed with the steady wind 
components and transformed to airplane axes for use in the calculations. 

The turbulence generated along each axis is shown in figure 2 and covers 
120 seconds of flight time. As the turbulence is the same from run to run, 
these curves show the turbulence the airplane is experiencing as it makes 
the landing maneuver. 


Acceptable Touchdown Conditions 

In the investigation reported herein, the landing was assumed to take 
place on a runway 3000 meters long and 50 meters wide. The desired touch- 
down point was on the center line 100 meters from the threshold of the 
runway. The values of 0 , > and Y> because of the reference 

used, indicate the errors from the ideal touchdown condition in which all 
these variables would have a value of zero. For the purpose of determining 
whether a landing was satisfactory, an arbitrary set of conditions was 
established. If a landing fitted within the following limits, it was said 
to be satisfactory: 

<j) = ±0.06 rad 
□ ±0.01 rad 

C-j = 0 to 500 m from touchdown point 

C 2 = ±10 m 

0 > y q >_ -1.98 x 10“^ rad 

0 < w -<1.0 m/sec 

z j _ 

0.017 ^ 9 ^ 0.061 rad 

The angle t gives the direction of the velocity vector when wind and tur- 
bulence are present or the direction of the airplane x-axis when there is 
no wind nor turbulence. It is always referred to as the heading angle. 


RESULTS AND DISCUSSION 
Autoland Longitudinal Control 

The longitudinal control of the basic system described in references 1 
and 2 is achieved by the following* three subsystems: 
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(1) The flight-path control system - this system receives guidance infor- 
mation from the guidance computer and controls the letdown of the airplane 

to the runway and, working with the sideslip control, provides turn coordination. 

(2) The automatic trim system - this system makes use of the stabilizer 
to trim continuously the pitching moment due to angle of attack. 

(3) The speed-control system - this system uses automatic throttles to 
control the speed of the airplane to specified values of approach speed and 
landing speed that are contained in the guidance computer and fed to the 
speed-control system as a function of flight-path angle and altitude. 

These three systems operated in combination to provide adequate automatic 
longitudinal control for the autoland system. Speed regulation was poor but 
sufficient for the landing maneuver. 

The results of further investigations on the flight-path and speed- 
control systems that improve the response of these systems and constrain 
the flight path to a specified plane are presented in the following sec- 
tions. The flight-path control system is discussed first. 

Flight-pat h control system .- The flight-path control system used in 
references 1 and 2 is essentially the same system. When a vertical dis- 
placement occurs as the result of wind or turbulence , the guidance does 
not attempt to return the airplane to the previous flight path, but it 
determines a new letdown strategy so that the airplane can fly to the land- 
ing point from its present position. 

Gain . change s and flap settings .- Further investigation showed that 
improved performance could be obtained by making minor modifications to 
the basic system. These modifications are to replace the term a with 
a sin and to change k^ from 1.0 to 0.8 in equation (17). In addition 
to these gain studies, it was found that lowering the upper limit on flap 
deflection helped the airplanes attain a proper pitch attitude for landing 
(ref. 1). In still air, limiting the deflection to 0.61087 radian improved 
the pitch altitude from -2.25 x 10“ 2 radian to 4.79 * 10~ 2 radian. Thus, 
limiting the flap deflection changed an unacceptable pitch attitude to an 
acceptable one. The maximum flap deflection was found to vary with wind- 
speed. For a wind of -25.8 m/sec (headwind) at 0.7584 radian to the runway 
center line the flap deflection used was 0.2618 radian so as to achieve 
a pitchup attitude at touchdown. As this wind vector was the only one used, 
the exact variation of the maximum flap setting with wind is not known. 

Letdown-system constraints .- The letdown system used in the basic auto- 
land system is constrained so that letdown does not start until the flight- 
path-angle command at the start of a landing maneuver is greater than 
-0.073 radian. The airplane maintains constant altitude flight until the 
letdown conditions are met. The ability of this system to hold altitude 
while making turns was demonstrated in reference 1 . The' introduction of 
winds and turbulence did not alter this ability. Flight-path-angle con- 
straints up to -0.10472 radian have been used successfully. 
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A flight-path constraint that forces the system to control the flight 
path to a specified plane was also investigated. Note that the flight-path 
constraint selected is only one of a number of possible constraints; for 
instance, a curve surface could have been used. The flight-path constraint 
used was selected because it was thought to be more generally applicble than 
the others. 

The plane used is defined by two lines that are perpendicular to the 
plane defined by the runway center line and the vertical to it. The first 
line passes through the present position of the airplane, and the second 
through the start-of-flare point. (See fig. 3.) The slope of the plane is 
given by 


t = tan-’ f - - 20 ) (18) 

\ R 11 ' ’” 5 ) 

The angle £ is continuously calculated until it is equal to a predeter- 
mined value of y called y„. When £ = y , letdown starts and y c is 
calculated by the formula 

Y c = Y p cos * (19) 

where 4 1 is the heading angle of the velocity vector with respect to the 
runway center line . The use of the heading angle in y c as a cosine func- 
tion is so that the airplane decreases its flight-path angle as it turns 
towards the perpendicular to the runway center line. 

In addition to equations (18) and (19), a logic block is necessary for 
those cases that have an initial heading angle greater than ir/2 radians. 

This logic holds y Q - 0 and prevents the calculation of £ (eq. (18)) 
until <0 and »P 1 ir/2. 

The altitude time history, altitude track, and ground track for a plane- 
constrained letdown are shown in figure 4. The ground track (fig. 4(c)) 
is the usual ground track for this initial condition, which is 
? 1 = -3000.0 meters, = -4000.0 meters, h = 200 meters, and = 0.0. 

Wind and turbulence were not used in this case. Two curves are shown for 
the altitude track and altitude time history. The curve with the plus 
marks is the desired flight path, and the unmarked curve is the actual 
flight path. The plus marks are also time marks and are spaced at 5-second 
intervals. As can be seen, the airplane when letting down is actually 
flying slightly above the desired flight path. This type of flight is 
acceptable as it is better to fly slightly above the glide slope than below 
it. In addition to the change in y c it was necessary to change the gain 
k^ in equation (17) from 1.0 to 2.1. The touchdown conditions for this 
case are 

♦ = -4.72 x io" 2 rad 
= 1 .895 x 10 -2 rad 
C 2 = -2.624 m 
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V. j, = 69.89 m/sec 


0 =3.09 x 10” 2 rad 
C-l = 2.83 x 10 2 m 
y = -8*761 x 10“^ rad 


= 0.061 m/sec 

j 

and all are within acceptable limits. Figure 5 shows the altitude time 
history and altitude track for plane-constrained flight at the same initial 
condition when wind and turbulence were added to the simulation. There was 
no discernible alteration in the ground track. Even in the heavy wind and 
turbulence used in this problem the system was able to control the flight 
path to within reasonable limits of the desired flight. The large excursion 
that starts at about 68 seconds is caused by the airplane entering a region 
of heavy horizontal turbulence and down blasts. (See fig. 2.) In spite of 
the severity of the turbulence the system was able to maintain control of 
the airplane and effect a reasonable landing. The touchdown conditions for 
this case are 


4) = 

1.391 X 

10 -2 

rad 

♦ = 

1.534 X 

10 -1 

rad 

C 2 = 

1.177 > 

< 10 1 

m 

11 

EH 

> 

7.421 > 

< 10 1 

m/sec 

0 = 

2.775 X 

10 -3 

rad 


^ b 9.873 m 
y = -1 .446 x 10“ 2 rad 

w z . = 1 .073 m/sec 
3 

Three touchdown conditions \p, and w z < are unacceptable. The cause of 


ij> and c ^ not meeting standards is probably the turbulence that hit the 
airplane at 60 seconds. As the landing took place 11 seconds after the 
onset of the extreme turbulence, little time was available to recover and 
correct path fluctuations caused by the turbulence. The condition w 

z . 


1.073 m/sec above the largest permissible value. This excess vertical speed 
is attributed to the down blasts and poor speed-control system. The speed 
was about 7.0 m/sec above the commanded value of 67 m/sec. If the commanded 


value had been met, 


w z # would have been below the maximum permissible 
3 


value for this parameter. 
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Speed-control system .- The details of the original speed-control system 
are given in appendix A. As indicated in appendix A, the gain has a 

value of 1.0. The curve labeled A in figure 6 shows the speed regulation • 
obtained with this system (see table I) in still air. The initial speed was 
77.12 m/sec and, as soon as the airplane started to let down, between 0.0 
and 1.0 second, a speed change to 72.43 m/sec was given. At about 40 sec- 
onds the speed entered the dead band of ±1.0 m/sec and remained within the 
dead band until a new speed command of 67.0 m/sec was introduced at about 
73 seconds. The control system was never able to bring the speed to within 
±1.0-m/sec dead band around the new command speed. The root-mean-square 
(rms) error in speed from 0 to 73 seconds was 3.12 and from 73 seconds on 
was 3.11. Curve B shows the effect of reducing the gain k^ to 0.45. 

The dead band was achieved at about 40 seconds but the large overshoot at 
the beginning followed by the undershoot caused the rms error to increase to 
3.52. However, better speed regulation was achieved. The dead band was 
achieved 2 seconds before touchdown, the end of the trace. In this case the 
rms error for regulation around 67 m/sec was 2.83. In spite of the poorer 
rms error over the first part of curve B, it was considered more acceptable 
than curve A because of the better speed control over the critical final 
part of the run. The landing conditions for the runs from which curves A 
and B were abstracted were acceptable. The speed response when wind and 
turbulence were added (see table I) is shown as curve C in figure 6. In 
this case, the speed never enters and holds within the ±1-m/sec dead band, 
and the time represented is shorter because the airplane landed short of the 
touchdown point. Because the excursions from the commanded velocity are 
less than in other cases at the beginning of the run, the rms error over the 
first 74 seconds is slightly smaller, 2.98; however, from that point to the 
end of the run the rms error is 7.80. 

Because of the poor response of jet engines to calls for thrust changes 
and because further changes In k^ give no significant improvement in 
speed response, it was decided to investigate the use of aerodynamic sur- 
faces as speed-control devices. 

Pure drag surfaces .- The first type of surface to be considered for 
speed control was a pure drag surface or speed brake. When this surface was 
used, the airplane was trimmed with the speed brakes undeflected. Deflec- 
tions of the speed brakes were used to control the speed of the airplane by 
changing the deflection angle of the speed-brake panels. It was assumed 
that the speed brakes were flat plates with a drag coefficient of 1.0 when 
they were deflected perpendicular to the plane of symmetry of the airplane. 
For intermediate settings it was assumed that the drag coefficients varied 
directly as the sine of the deflection angle. Details of the system used to 
control these speed brakes are given in appendix B. The ability of the drag 
flaps to control the speed of the airplane is compared with that of the auto- 
matic throttle in figure 7. The automatic-throttle speed control is repre- 
sented by curve C and the drag flaps by curve D. Over the first 74 seconds 
the rms error for the drag flap was 1.12 as compared with 2.96 for the auto- 
matic throttle. During the time from 74 seconds to touchdown, the rms 
errors were 2.42 and 7.80, respectively. The touchdown conditions for drag 
flaps were acceptable except for the vertical touchdown speed, which was 
marginal at 1.39 m/sec. This higher than acceptable touchdown speed proba- 
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bly results from the high turbulence that occurred from 68 to 77 seconds. 

A comparison of curves C and D of figure 7 shows that the drag flaps provide 
faster and more precise control of the speed than do the automatic throt- 
tles. However, in order to accomplish the speed regulation, the drag flaps 
had to provide, based on wing area, an increase of 0.11 in the drag coeffi- 
cient. Under the assumptions concerning the drag flaps, this means that the 
drag flaps would have to have an area equivalent to the area of the wing 
flaps. Because of the size equivalence it was decided to investigate the 
use of wing flaps as a speed-control device. 

Wing flaps .- Except for the change in constants (see appendix B), the 
control system used for the wing flaps was the same as that used for the drag 
flaps. The wing flaps were at an initial deflection of tt/ 12 radians and 
the thrust trimmed . This thrust was then held constant throughout the land- 
ing maneuver. The speed regulation obtained from the use of wing flaps is 
compared with that obtained from automatic throttles in figure 7- Curve E is 
for the wing flaps and curve C is for the automatic throttle. The rms over 
the first 74 seconds was 0.51 for wing flaps and 2.96 for the autothrottle. 
Over the critical final part of the landing maneuver, the rms errors were 
1.58 and 7.80, respectively, for the wing flaps and automatic throttles. The 
rms error in speed obtained for the wing flaps represents an improvement over 
that obtained for the speed brakes. The wing-flap speed-control touchdown 
conditions given in table I are within acceptable limits. An examination of 
the accelerations along the x body axis of the airplane showed that the 
maximum acceleration was -0.22 g unit (1g □ 9-80665 m/sec ) and occurred just 
after the speed change commands were given. This acceleration is a little 
high for passenger comfort but could be reduced by restricting the flap- 
deflection rate. The most significant difference is the time taken by the 
speed brakes to reduce the speed after a command was given. The speed- 
reduction commands occurred at points 0.001 second and at 74 seconds. After 
the first command, it took the speed brakes 14 seconds to achieve a speed 
level that the wing flaps achieved in 4 seconds. After the second speed com- 
mand, the speed brakes were not able to reduce the speed to 67 m/sec before 
touchdown, whereas the wing flaps were able to reduce the speed to 67 m/sec 
in 6 seconds. The delay in speed regulation associated with the speed 
brakes at the start of the landing maneuver is tolerable, as it occurs far 
from touchdown and between 14 and 74 seconds. Differences in speed regula- 
tion are negligible . However, the delay in achieving 67 m/sec is more crit- 
ical, as it occurs during the last few seconds of flight. 

The motions of the speed brakes and wing flaps, when used to control 
speed, are about the same; both are responding to turbulence. See figures 8 
and 9 for time histories of the speed-brake and wing-flap motions. The most 
erratic motion occurs around 70 seconds when the turbulence is most severe. 
Although the motions are equally erratic, the damping gain required with 
speed brakes was about 2.5 times that required for wing flaps. 

The results obtained from the three methods of speed control, autothrot- 
tle, speed brakes, and wing flaps, indicate that wing flaps, because of the 
precise speed control, constitute the most desirable method. 
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Wing-flap speed control was tried for other wind conditions than those 
given in this section. The speed control was equally good under all condi- 
tions investigated, but it was found that flap-deflection angles had to be 
increased so as to obtain good speed regulation from a windspeed of 
25.80 m/sec to zero windspeed. 

Effect of speed__control __on touchdown__condit ions . - An examination of the 
touchdown conditions for the various methods of speed control given in 
table I shows that case E has the best touchdown conditions of the group. 

The ground tracks for cases C and E are shown in figure 10. The only dif- 
ference between the systems that produced the ground tracks is the method of 
speed control. When wing-flap speed control is used, the first turn is 
tighter and this influences the ground track and letdown. Note that curve E 
alines with the runway sooner than does curve C, and there is a short period 
of straight flight before touchdown. The effect cf the tighter first turn, 
caused by better speed control, causes the small differences in the subse- 
quent ground track. This effect is apparently the main factor in producing 
the better touchdown conditions of case E. 


Lateral Guidance Considerations 


At the beginning of this paper material was presented that modified the 
lateral guidance so that initial heading angles greater than tt/2 could be 
used. An initial condition ^ = -1012 meters, ^ n -4000 meters, 
h = 540 meters, and i|>. = it was selected to demonstrate the adequacy of 
these changes. As wind and turbulence were used for this check, the heading 
angle of the velocity vector was less than ir radians. As can be seen from 
figure 11, both the ground track and altitude were good. The landing condi- 
tions for this case were not satisfactory, the overshoot was too great, and 
both and 6 were outside of the limits. The poor landing conditions 

were due to the shallowness of the glide slope, which was about 
-0.01745 radian at the start of the flare. When the glide slope was 
increased so that at the start of the flare it was -0.05044 radian, all 
landing conditions were satisfactory. The glide slope may be increased by 
increasing the initial altitude or by restricting the letdown until the 
heading angle is equal to or less than tt/2. The former method was used in 
this study. 


In the landing maneuver shown in figure 11, the turn is simple, that is, 
the roll angle does not change sign; therefore for the guidance algo- 

rithm (eq. (16)) is determined by the expression 


k 


1 = 


»c - *0 - Mo 
6 Ac " e AA 


This expression is used for the second calculation of in appendix B of 

reference 2. 


Avoidance of restricted areas .- In the context of this paper, restricted 
areas are interpreted as mountains, water towers, buildings, and ground 
areas over which flight is not permitted. The basic system reported in ref- 
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erences 1 and 2 and used in this paper was to be used over the last 
120 seconds of flight, that is, the final approach and touchdown. During 
this period it was assumed that all restricted areas would have been passed. 
However, restricted areas could be handled by assuming that the desired 
heading angles are supplied from ground control and by using the same turn 
algorithm (eq. (16)) as normally used for landing approaches. In this case, 
k. is set equal to zero and kp is positive. The values of k^ and kp 
would be maintained until the final approach was started. At this time k^ 
would be determined by the method given in appendix B of reference 2 and the 
sign of k 2 changed so that kp is now negative . 

Landing maneuver sensitivity to gain k^ .- When the initial position of 

the airplane is such that < |cp| the execution of the landing maneu- 

ver is very sensitive to the value of k^ in equation ( 1 6 ) . The sensitiv- 
ity increases as the wind decreases and is maximum at the zero wind condi- 
tion. Table II gives the landing conditions for the initial condition 
g. = £p = -4000 meters, h = 540 meters, and ^ = 0.0 for various values 
of k.fl). The value of k^l) was specified and the second value was cal- 
culated in the usual manner. No wind or turbulence was used in these calcu- 
lations. In the following discussion two values of k 1 will be referred to 
k^d), the k^ calculated at the start of the landing maneuver, and k 1 (2), 
the value of k 1 calculated at the completion of the first, that is, as <f 
passes through zero. The points of calculations of these gains are shown on 
a typical S-curve ground track in figure 12. Inspection of these data 
(table II) shows that a value k^l) between 1.55 and 1.60 gives the best 
all around landing, whereas for k^l) = 1.75, a value determined by the 
guidance algorithm, the heading angle is outside acceptable limits. When 
k. was increased to 1.80, the heading angle was larger and ? 2 increased 
from -2.15 meters to -6.08 meters. In neither case has the airplane reached 
the runway center line; whereas in case G, k^(l) a 1.55, while not reaching 
the center line was almost centered on the runway. The ground tracks and 
heading-angle time histories for cases G and K are shown in figure 13. The 
small value of k^(i) in case G increases the length of time during the 
first turn so that the airplane velocity vector is closer to the perpendicu- 
lar to the runway center line, and when k^(2) is calculated, it has a 
smaller value and lengthens the time devoted to the second turn. The total 
increase in time for case G over case K was 5.7 seconds, enough to give bet- 
ter landing maneuver performance. One way to lengthen the time to land and 
consequently the length of the ground track without changing the calculated 
value of k^(1) is to add the integral of the error to the turn command. 
When this is done, equation (16) becomes 

^c = k 2 [j^c “ ♦ - + k i ( e AC " 0 AA)] 



where k-^ is the integrator gain. Case M of table III gives k^(1), 
kp(2), and the landing conditions obtained where an integrator was used. 
Case M is the same as case K except for the integrator. The gain k IG was 
set to 0.75 during the interval between points A and B on the ground track 
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shown in figure 12 and to 0.0 from point B to touchdown. The ground track 
and heading-angle time history for case M are shown in figure 13. The 
curves for case M are almost the same as those for case G, as are the land- 
ing conditions given in tables II and III. When a wind of -25.8 m/sec at 
tt/ 4 radians to the runway center line was included, no integrator was neces- 
sary (see case N in table III); however, when the wind was reduced to 
-12.90 m/sec, case 0 in table III, an integrator was again necessary. How- 
ever, only half the amount of gain for the still-air case was required. 

The value of ( 1 ) was about halfway between those for cases G and M. The 

ground tracks and heading-angle time histories for cases N and 0 were prac- 
tically the same as those for case M. From these results it is apparent 
that the value of k^ is dependent on the value of k^ , which is a func- 
tion of airplane position and windspeed at the time k^ is determined, thus 
making k^ G dependent on those quantities. When k^ G was determined by 
the expression 

k lG = 3.755k., - 5.819 (21) 

satisfactory results were obtained. The gain k determined by equa- 
tion (21) was subject to two limitations: The first is for kj G < 0, then 

kjQ = 0.0; and the second is, if k^(2) has been calculated, k = 0.0. 

The first condition prevents a value of k^ from being used where an inte- 
grator is not needed as in case 0 of table III. The second condition 
restricts the use of the integrator to the first turn of the landing trajec- 
tory, for example between the points labeled A and B on the ground track of 
figure 12. Eauation (21) was evaluated over a wide range of initial condi- 
tions and was found to give satisfactory results. 

Eauation (21) applies to cases where automatic throttles were used for 
speed control. A similar problem exists when flaps are used as a speed- 
control device; however, different constants would have to be used in 
equation ( 21 ) . 


CONCLUDING REMARKS 

The letdown and speed-control subsystem of the longitudinal control sys- 
tem of the possible automatic landing system described in NASA TN D— 7 6 1 1 and 
NASA TN D— 797 1 have been studied. In the case of the letdown system a 
method for constraining the airplane to a letdown in a specified plane is 
presented. Constraints that control the steepness of the descent path of 
the original letdown system are also given. 

The effectiveness of aerodynamic surfaces instead of automatic throttles 
as speed-control devices was investigated. In general, drag flaps (speed 
brakes) or wing flaps provided more precise speed control than the automatic 
throttles. When the aerodynamic surfaces were used, the vertical touchdown 
speed was higher than with the automatic throttle , but within the design 
limits of transport aircraft. Because of less flap activity, wing flaps 
seem to be preferable to drag flaps. 



In addition to the work on the longitudinal control system, modifica- 
tions to the lateral control system to permit zr-radian turns are presented. 

A method of using integrators to reduce the sensitivity of the lateral guid- 
ance system to the gain k^ (variable gain for the pseudo-radar azimuth- 
angle error of the turn-control algorithm) in the roll command was 
developed . 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 17 , 1976 
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APPENDIX A 


AUTOMATIC-THROTTLE SPEED-CONTROL SYSTEM 

Speed control was achieved through the use of a digital thrust control 
system. The following eouation was used to compute change in thrust AT 
required to eliminate a speed error 

AT = k 6 ( V c ~ V t) " k i5 m ^ (A1) 

where 


and 


k 6 = V c (28.131 + 58.9736 f ) (A2) 

If the aT calculated by equation (A1) is added to the present thrust 
T, the new total thrust is T + AT. The thrusts T and T + AT were 
looked up in table I of reference 1 , which gave information on how the 
thrust changed with time. The output of this process was At, the time it 
would take the engine to change its thrust from T to T + AT. The change 
in thrust was added linearly over the time period At, so when At seconds 
had passed, the engine thrust was T + aT. After the thrust was changed, 
the thrust controls were shut down for 4 seconds. At the end of 4 seconds, 
V,p and u were sampled to determine whether they were within the control- 
system dead band. The dead band for V T was ±0.1 m/sec and for u was 
±0.1 m/sec . If the errors were within the dead band, no more changes were 
made. However, if one or both errors were outside the dead band, thrust 
changes were made until the error entered the dead bands. 
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II 


APPENDIX B 

FLAP SPEED-CONTROL SYSTEM 

The flap control system used for speed control is shown in the following 
block diagram: 



where k^g = 1.461 and = -0.03 when wi ng flaps are „sed, k^g = 1.190 

and k^ □ -0,075 when drag flaps are used, ar i t - C.C^'L for both 
cases . 

Flap deflection is positive from the zerc f ] \o 00s it . Fecause 
represents the estimated change in flap deflection to ^iv. , ‘ lr ' desired speed 
change, the current position of the servo mus" t-~ ^dde^ r , to obtain 

. The state-variable convolution method w*s ur-ed tc r f \ > -| she inte^ra- 
1 c 

tions in this control system. 
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TABLE I.- LANDING CONDITIONS FOR SPEED-CONTROL CASES 


[Initial condition: ^ = -3000.0 m, Cg = -4000.0 m, h = 540.0 m, □ 0.0; wind, 

-25.80 m/sec at ""74 rad to runway center line; patchy turbulence (see fig. 2)] 


Case 

Variable^. 

A 

(a) 

B 

(b) 

C 

(c) 

D 

(d) 

E 

(e) 

<t> , rad 

-6.91 x 10' 2 

-6.27 x 10 -2 

-1.36 x io" 2 

1.52 x io' 2 

3.23 x 10" 3 

4 1 , rad 

2.08 x io -2 

1.932 x io -2 

1.92 x io" 2 

1.32 x io -2 

-1.52 x IQ -3 

^2 » m 

2.29 

1.83 

-7.13 

-8.93 x 10' 2 

9.84 x 10" 2 

c 1 , m 

4.83 x 10 2 

4.79 x io 2 

-9.544 x io 1 

3.34 x 10 2 

2.894 x 10 2 

y, rad 

-7.71 x 10" 3 

-7.74 x 10" 3 

- 3.60 x io" 2 

-2.065 x 10" 2 

-6.58 x 10" 3 

9 , rad 

3.55 x io -2 

4.79 x io" 2 

-7.34 x io" 2 

2.75 x io' 2 

1.66 x io" 2 

w z , m/sec 

5.32 x io -1 

5.20 x io -1 

2.64 

1.39 

.44 

J 

V T , ra/sec 

6.895 x io 1 

6.713 x io 1 

7.332 x io 1 

6.736 x IQ 1 

6.697 x 10 1 

t , sec 

83.39 

84.06 

1 

78.63 

86.65 

86.45 


a Autothrottle speed control; acceleration feedback erain, 1.0; no wind; no 


turbulence . 

D Autothrottle speed control; acceleration feedback gain, 0.45; no wind; no 
turbulence . 

c Autothrottle speed control; acceleration feedback gain, 0.45; wind and 
turbulence . 

a Drag-flap speed control; acceleration feedback gain, 0.075; wind and turbu- 
lence; no thrust cut. 

e Wing-flap speed control; acceleration feedback gain, 0.01; wind and turbu- 
lence; no thrust cut. 



TABLE II.- EFFECT OF 

[Initial condition: ^ ° -^000.0 m, h 


Case 

Variable. 

F 

G 

H 

^(D 

1.50 

1 .55 

1 .60 

k 1 ( 2 ) 

1.472 

1.52 

1.57 

<t> , rad 

2.57 x 10" 2 

1.87 x 10" 2 

9.34 x 10“ 3 

0 , rad 

2.59 x 10“ 2 

2.42 x 10 -2 

1.76 x 10 -2 

, rad 

m 

o 

X 

Lf> 

-3.53 x 10“ 3 

-3-79 x 10" 3 


154.2 

126.1 

141 .5 

C 2 > ® 

-3.22 

-.842 

.474 

Y , rad 

-7.14 x 10“ 3 

-5.75 x id" 3 

-5.52 x 10" 3 


TOUCHDOWN CONDITIONS 


5^0 m, = 0.0; no wind or turbulence] 


I 

J 

K 

1 .65 

1.70 

1.75 

1 .62 

1.67 

1.72 

3.09 x 10 -3 

-2.53 x 10" 3 

-4.41 x 10“‘ 

o 

X 

o 

1 

no 

4.42 x 10 -2 

4.75 x 10 -: 

6.37 x lo" 4 

8.97 x 10" 3 

2.26 x lo - ' 

>2.4 

53.7 

50.1 

.724 

.444 

-2.15 

on 

l 

o 

X 

o 

C\J 

-1.53 x 10" 3 

-6.83 x 10“ 



TABLE III.- EFFECT OF INTEGRATORS IN TURN ALGORITHM 


[Initial condition: c 1 = c 2 = -4000.0 m, h □ 540 m, ^ = 0.0] 


Case 

M 

N 

0 

Variable 




Integrator G 

0.75 

0.0 

0.375 

k-l ( 1 ) 

1.75 

1.58 

1.65 

k-| ( 2) 

1.55 

1.55 

1.57 

4> , rad 

1.297 x 10- 2 

2.33 x 10"2 

1.729 x 10-2 

e , rad 

2.02 x 10“2 

-4.59 x 10-2 

7.34 x 10~3 

\|j , rad 

-3.9 x 10"3 

3.038 x 10"3 

2.741 x 10 - 3 

C 1 ? m 

138.7 

74.75 

52.9 


-.077 

-.493 

-.5 

Y , rad 

-5.89 x 10-3 

-1.26 x 10"2 

00 

1 

o 

X 

[>- 

GO 

O 

1 

V T , m/sec 

70.14 

71.94 

70.22 

, m/sec 

v 

.413 

.906 

.763 

t , sec 

94.5 

92.2 

92.7 

Wind, m/sec 

0 

-25.80 

-12.90 

Turbulence, m/sec 

0 

Patchy 

0 
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Figure 2 
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(a) Turbulence along X-axis. 
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(b) Turbulence along Y-axis. 
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(c) Turbulence along Z-axis. 

Turbulence used in study, referred to as patchy turbulence, 
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(b) Altitude track. 

Figure 4.- Altitude time history, altitude track, and ground track for 
plane-constrained case; no wind or turbulence. Curve with plus marks 
is desired flight path; unmarked curve is actual flight path. Plus 
marks are spaced at 5-sec intervals. 
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Figure 4.- Concluded. 
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(a) Altitude time history. 
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(b) Altitude track. 



Figure 5.- Altitude time history and altitude track for plane-constrained 

case with wind and turbulence. 
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Figure 6.- Time history of speed for automatic- throttle curves A and B; no wind and no turbulence. 

Curve C; with wind and turbulence. 
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Figure 7.- Comparison of speed control in wind and turbulence for curve C (automatic throttles), curve D 

(drag flaps), and curve E (wing flaps). 






5 r m 

Figure 10.- Effect of speed-control method on ground track in wind and 
turbulence. Curve C (autothrottle) and curve E (wing flaps). 
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Figure 11.- Concluded 
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